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1. INTRODUCTION 

Nowadays, the increasingly high number of mobile users has put pressure on current 
communication systems adopting such techniques as time-division multiple access (TDMA) and frequency- 
division multiple access (FDMA). The overload is especially intensified in the new 5th Generation system. In 
order to tackle the problem of massive users and increase access capacity, beam division multiple access 
(BDMA) has been proposed to the system [1], [2]. The BDMA scheme basically consists of three steps: i) 
base station (BS) acquires channel coupling matrices (CCMs) from all user equipments (Ues) in its cell, ii) 
user selection based on CCMs, then beams at BS are group into separated sub-sets, and iii) transmitting pilot 
and data in the uplink and downlink. A selecting scheme within non-overlapping beams is proposed to 
BDMA based on decomposition of the MU-MIMO channels into multiple single-user MIMO channels [1]. 
As a result, the overhead of channel estimation, as well as the processing complexity are improved 
significantly at transceivers. Meanwhile, Jiang et al. [2] apply lyapunov-drift optimization framework for joint user 
scheduling and beam selection method, thereby obtaining an optimal scheduling policy in a closed form. 

The average channel capacity of wireless communication systems is heavily affected by fading 
environments [3]-[6]. Especially, channel capacity degrades significantly it’sin the fading correlations 
enivroments [7]. Massive multiple-input multiple-output (MIMO) techniques, for example, hybrid 
beamforming architecture, were developed to attain the improving channel capacity of system in rich 
scattering environments. Hybrid beamforming was deployed widely for millimeter wave (mmWave) massive 
MIMO transmissions [8]-[11]. This architecture is considered the most effective solution to power 
consumption and production costs problems of fully-digital architecture based on an additional layer of 
analog phase shifters [12]. Moreover, the hybrid beamforming also achieves efficiency approximating to the 
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performance of digital solutions but at a much lower cost [13]. Transmit beamforming strategies for the 
cellular downlink become the attractive topic [14]-[19]. There are several previous works proposed such 
efficient beamforming strategies as zero-forcing beamforming (ZFBF) [14], [20]-[22], orthogonal 
beamforming (OBF) [23], [24], and orthogonal random beamforming (ORBF) [4]. If the number of mobile 
stations is higher than the number of antenna, user selection technique is essential to selects a group of users 
for simultaneous transmission [25]. Works in [20], [21] proposed two selection algorithms for MIMO which 
are semi-orthogonal user selection (SUS) and greedy user selection (GUS). In [21], by combining ZFBF with 
SUS (ZFBF-SUS), the optimal sum rate could be achieved asymptotically with fully-digital, as the number of 
users approaches infinity. 

In this work, we study a low-complexity downlink-beamforming algorithm for step 2 in BDMA 
scheme such that maximize sum capacity for the practical case when the number of downlink users exceeds 
the number of antennas. Therefore, we propose a multiple access scheme that employs BMDA design with 
hybrid beamforming structure and simple user-beam schedule, which not only reaches asymptotically 
optimal in the sum-rate, but the average multiuser interference (MUI) energy can be minimized. In each time 
slot, before the antenna weight vectors are updated base on ZFBF, the base station will select those mobile 
stations such that the interference between them is minimal. To be specific, the system aims at 
opportunistically schedule users whose steering vectors are as orthogonal as possible to each other. 
Subsequently, this scheduling scheme will reduce the average multiuser interference MUI energy over all the 
mobile stations (MSs). 

The content of paper is organized as 5 sections. In section 2, we present the problem formulation, 
and ZF with user selection algorithm is considered in section 3. Results and discussions are given in section 4 
to demonstrate the effectiveness of our algorithm. Some conclusions are given in the last section. Notation: 
uppercase boldface letters and lowercase boldface are described correspondingly for matrices and vectors. 
The operators||. ||, (.)7, (.)” and |. | denote norm, transpose, Hermitian conjugate and modulus, respectively. 


2. PROBLEM FORMULATION 

Generally, the hybrid beamforming architecture could be divided into fully connected and partially 
connected one [26]. While each radio-frequency (RF) chain of a fully connected architecture is mapped with 
all antenna elements by phase shifter network, a partially connected hybrid transmitter connects each RF 
chain to an antenna sub-array. In this paper, we conducted our proposed algorithm with fully connected 
hybrid architecture as shown in Figure 1. 
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Figure 1. A fully connected hybrid beamformer 


We consider a massive mmWave MIMO system with Nr transmitted antennas and Ngp RF chains at 
BS serving Ny mobile station (MS) with Ng received antenna at each one. In this work, we concentrate on a 
single block transmitted to Ny users. The resulting signal x of dimension Ny X 1 after hybrid beamforming 
system is modeled as in (1): 


X=FppXFep XS=FXS (1) 


where FRf = [f Rr, PRE > fupe RF] is analog beamforming matrix with fu prp E CNTX1 being the u-th analog 
beamforming vector, the digital beamforming matrix Fgg = [f sB fo BBs) fer BB| with fugg E CNs*? is the 
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digital beamforming vector for the u-th user, and sy is the block data of u-th user. The received signal of the 
u-th user can be expressed as (2): 


Yu = Hy: fu,RF : fue ‘Sy + > H,, ' fu RF : fi eB Si e n, 7 
UE 
desired itu noise ( ) 


interferece 


where y, E CNR*1, H, E CNR*NT is the MIMO channel matrix between the transmitter and the u-th receiver, 
and ną E€ CR*1 is the complex additive white gaussian noise with zero mean and variance of g? for u-th 
user. At u-th MS, the beamforming weight vector w combined the RF wer and baseband wg,z precoding 
vectors. The decoded signal of the u-th receiver is given by (3): 


~ aa H H ~ 
Su ai Wy Hufu rr fu ee su F Wy Ny (3) 


where ñ, is the sum of complex additive white Gaussian noise (AWGN) and interference from other users. 
The u-th user’s channel model of mmWave systems can be modeled as in (4) [27], [28]: 


N,N, 
J aS ae ar(Pip Pu) a 4 (Pip Pus) (4) 


where N, is the number of scatters of the u-th user’s channel, a, is the complex path gain, arC; p Puy) and 


ar(o. p ONE are received and transmitted steering vectors. In this work, a N} X Nw uniform planar array 
(UPA) is used, the steering vector can be rewritten [29]: 


alọ, p) = a [1, eskd(sin ¢ sin p+cos p) — eskd((N; —1) sino sin p+(Nw-1) cos oy (5) 
T 


2n . ; ) i 
where k = = is the wavenumber and d is the distance between two adjacent element antennas. Then, the 


signal-to-interference plus noise ratio (SINR) of user u" is (6), 


|Ha ° furr ' fusa) 


SINR, = —— (6) 
Disul Hu ' fu RF Gal +0? 
The beamforming cost function can now be formulated as (7): 
Nu 
par 2 log2(1+ SINR,,) subject to ||Frp. Fggll? = (7) 
RFF BB 


3. PROPOSED ZF WITH GREEDY ALGORITHM 
3.1. Proposed algorithm 

We assumed that BS has the perfect channel state information (CSI). Therefore, angle of arrival 
(AoA) and angle of departure (AoD) information {Ø}, v/,, pt, p£} is perfectly known. As a result, optimal 
analog beamforming vector of u-th user could be computed (8): 


furr = ar(ġi Py), u= 1+Ny (8) 


The optimization problem in (7) is rewritten (8): 


Nu 
max) log2(1 + SINR,,) subject to||FRr-Feall* = 1 (9) 
BBs 
From (3), set g? = wi! x H, x fire» then G = ese E TR According to a zero-forcing approach, in 


case of Ny > Ny, the optimal DBF matrix could be calculated as F pg = Gt = GĦ (GGĦ)~1. Normalized 
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; . F 

baseband precoder F gz is obtained Faz = EET For Nr < Ny, to select subset users, we proposed a low- 
RFF BB 

complexity algorithm, named ZF with user selection (ZFUS), as summarized in next. 


Algorithm 1. ZF with user selection (ZFUS) 

Initialization: N,Ny, Nr, Nr, Neg, SNR, AoA and AoD {p} o5 bt, pt} 

Beam search 

— Find the user with the largest channel gain among Ny users: 
o sı =argmax H,,Hy 


-— Set its steering vector as first element of signal subspace: A= [a,(@f,9%)] 
- n=; 

Select users 

= while(n<N—1) 


n=nt+1 
The projection space: Py = AAt 
fori =1:Ny 
P(i) = [IPa * ar(t ÐI 
end 


Find i = argmin{P(i)} 
Update signal subspace:A = [ar(o, @5),..., ar (of, o] 
end 
Obtain list of Nselected users 
Analog beamforming: Fer = [ar(ot,9),..., Ar (Phar Prize) | 
Digital beamforming: Fgg =Gt=G4(GG4), where G = [g,,g>,..., gylt 
End 


There are two main steps in user selection: beam search and select user. In the first phase, based on 
broadcast beam at BS side, we can find out UE with largest channel. Then, we set its steering vector 
ar({, p£) as first element of signal subspace A. For the next phase, by comparing orthogonality between 
user’s steering vector and projection space, the algorithm will select (N — 1) users among the rest of 
(Ny — 1) users. 


3.2. Computational complexity analysis 

In this subsection, the computational complexity of the ZFUS algorithm is analyzed. The 
implementation of ZFUS consist of two phases: User selection and beamforming vector calculation. Because 
the latter requires only N x N matrix inversion to calculate beamforming weight, we focus on the complexity 
of user selection phase. In the beam search step, we need one 2N; matrix-matrix multiplication per user. 
With Ny users, we need 2N;Ny matrix-matrix multiplication. The complexity of the selecting user step can 
be expressed as )it_,k (8Nrk? — Nrk +2(Ny — 1)Nrk) = X} k (8Nrk? + (2Ny — 3)Nrk). As a 
result, the number of computations of running the ZFUS algorithm is 2N;Ny + YRIi(8Nrk? + (2Ny — 
3)Nrk). 


4. SIMULATION EVALUATION 

The performance of proposed ZFUS algorithms will be illustrated numerically in terms of the 
average MUI energy over all the MSs, system sum-rate capacity conditions used in our simulation are: i) the 
inner distance of two adjacent antennas is 0.54, ii) the number of total users Ny = 20 users with a 2x2 UPA 
(Ng = 4) and the number of selected users N = 4, and iii) the number of realizations is 500 for each 
experiment. Figure 2 shows the system sum-rate capacity versus SNR in 256 X 16 system with uniform 
planar array at both the transmitted and received side. It can be seen that the performance of ZFUS is 
approached to fully digital precoding system with a small gap. When compared to analog beamforming 
(ABP), Figure 2 shows that the ZFUS significantly improves the performance of system. 

Next, we investigated the performance of three architectures following the number of base station 
antennas. Figure 3 shows that spectral efficiency (SE) values of three systems rise rapidly as the number of 
antennas increases from 16 to 250 elements. Though the SE of ZFUS is always higher than ABF about 4 
bps/Hz, it could not reach asymptotically the fully digital system. 

Finally, we study the Average of MUI energy over all the UTs versus number of BS antennas in 
Figure 4. It’s clear that ZFUS can reduce MUI more effectively by selecting the user group with orthogonal 
steering vector to each other. While the MUI value of the ABF is always greater than 0 dB, the ZFUS 
maintains the MUI value below -15 dB. As the number of antennas increases, both the MUI value of ZFUS 
and ABF tend to decrease. 
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Figure 2. Comparison of sum-rate capacity versus SNR with N_T = 16 x 16 
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Figure 3. Comparison of sum-rate capacity versus number of BS antennas 
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Figure 4. Average of multiuser interference at BS 
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Compared to several relevant papers in Table 1, the performance of our proposed is better under the 
same architecture and SNR. To reduce the computational complexity, Singh and Ramakrishna [30] use 
subconnected structure and codebook method, this leads to decrease the spectral efficiency of system. In 
contrast, though works in [31], [32] deployed the fully structure and zero-forcing method same as our work, 
their spectral efficiency values are still lower than our work. It proves that the proposed user schedule help to 
improve the performance of system. 


Table 1. The comparison of our proposed and previous works 


Paper Architectures Antenna Array Algorithms SE (b/s/Hz) 
N =2x8 
E ABF: codebook-based 
[30] Sub-connected Ng N, 2 A : DBF: a specified codebook of matrices 10@SNR 10 dB 
= 
MES Oe ABF: Hermitian of downlink channel matrix 
[31] Fully connected Np = 1 i À 1.37@SNR 0 dB 
DBF: ZF 
Ny = 2 
Nr = 64 i 
e ABF: beam steering codebook 
[32] Fully connected ae = T DBF: ZF based on perfect CSI 6.8@SNR 0 dB 
= 
Nr = 16 x 16 nes 
My Work Fully connected Np = 4 ant eee vector 10@SNR 0 dB 
N =4 DBF: ZF 
Woe 


5. CONCLUSION 

A new user selection algorithm is proposed for BDMA scheme of massive MIMO system in this 
work. Based on deploying Zero-Forcing and new schedule in hybrid beamforming, the multiuser interference 
(MUI) is significantly reduced when we increase the number of BS antenna. Simulated results confirmed that 
proposed solution can attain spectrum efficiency (SE) as good as only digital beamforming. 
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